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Abstract 



The time evolution from coherent states to squeezed states of high density 
excitons is studied theoretically based on the boson formalism and within the 
Random Phase Approximation. Both the mutual interaction between excitons 
and the anharmonic exciton-photon interaction due to phase-space filling of ex- 
^r) • citons are included in consideration. It is shown that the exciton squeezing 

depends strongly on the exciton density in semiconductors and becomes smaller 
with increasing the latter. 
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Recently the quantum optical properties of the quasiparticles in condensed mat- 
ters have been investigated intensively [1-10]. In particular, the nonclassical properties 
of exciton-polariton were considered in [1-3], phonon-polariton in [4], and squeezed 
states of excitons were considered in [5-7]. Of particular interest are also the optical 
q ! nonlinearities emerged from high external excitation which affect almost all the optical 

properties of semiconductors. At high excitation a great number of excitons can be pro- 
duced and there arise some kinds of anharmonicities such as nonlinear exciton-exciton 
interaction and the anharmonic exciton-photon interaction due to phase-space filling of 
excitons. An attempt to consider the effect of these anharmonicities on squeezed states 
of excitons was first made in [5-7]. However, in [5,6] to solve the problem the authors 
had to resort to rough approximation neglecting many terms in Hamiltonian which ac- 
tually should be retained. It is well known that for the high density systems where the 
fluctuations of the occupation number can be neglected, the Random Phase Approxi- 
mation (RPA) [11] is the most suitable, and one can safely apply it to study nonlinear 
optical properties of the system. In this paper, as in [7], based on the boson formalism 
and within the framework of RPA, we study the effect of both above-mentioned two 
kinds of anharmonicities on the squeezing properties of high density excitons. It will 
be shown that the squeezing degree of excitons depends strongly on the exciton density 
in semiconductors, and moreover, it decreases with increasing the latter. 
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We shall start with the following effective boson Hamiltonian describing the exciton- 
photon system at high excitation [9] 



H eff = E ex a + a + uo^c+c + fa + a + aa + [g Q ca + (l - \a + a) + h.c], (f ) 

where a + (a),c + (c) are the creation (annihilation) operators for exciton and photon; 
E ex , and uj 1 are the energies of the exciton and photon, respectively, g is photon- 
exciton coupling constant. The third term on the right-hand side of (1) describes the 
repulsive exciton-exciton interaction, and / is given by 

/ = ~ e v(k - k 1 ) [\m\ 2 Mk')\ 2 - \m\ 2 \m^k')\] . (2) 

V k,k' 

The fourth term on the right-hand side of (1) represents the anharmonic exciton-photon 
interaction due to phase-space filling of excitons, and A is given by 

a = ^EM*)IV(*)/Ev(*)- ( 3 ) 
v k k 

Here V is the quatization volume, f(k) - is the Fourier transform of the exciton envelope 
function for electron-hole relative motion, V(k — k') - is Coulomb potential. 

We shall not detail here the derivation of Hamiltonian (1) reminding only that it is 
obtained by operating Usui transformation [12] on the Hamiltonian of the interacting 
photon-electron-hole system. This procedure was originally developed by the authors 
of [13,14] for boson formalism of high density exciton theory. 

Let us consider the following Green's functions 

Gi(M) = 0(t) < [a(t),a + (o)\ >=« a(t),a + (o) », (4) 

G 2 (k,t) = «c(t),a + (o) » . (5) 

The equations of motion for the Green's functions (4,5) can be derived in standard 
way. Then using the RPA-type factorization which is reasonable for high densities, we 
obtain the closed system of equations, the solutions of which give us the spectrum of 
the exciton-polariton, as follows 

e 1>2 = 1 -{E* + ^) ± y(E*-^) 2 + 4g* 2 , (6) 

where 

E* = E ex + 2n ex Vf, (7) 

g* 2 = g 2 (l-n ex V\), (8) 
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and n ex is the density of excitons created by laser field in semiconductors. 

The polariton whose spectrum is (6) can be now expressed in terms of photon c 
and "real" exciton a operators, as follows 

B,i = UiC + v^, (9) 

where 

Ul = [l+g*\ ei -E*)- 2 }- 1/2 , (10) 

Vi = <?*|q-£*|-V (11) 

To analyze the time evolution of the squeezing properties of the exciton field we 
introduce two Hermitian quadrature exciton operators 

X(t) = ±[a + (t)+a(t)], (12) 

P(t) = *-[a + (t)-~a(t)}- (13) 

From (9) and the exact time- dependent expression for operator £>«(£) one can obtain 
the following equation 

a(t) = E [uiVic(o) + vta{o)] e~ ie ^ . (14) 

i=l,2 

Since the squeezed states of exciton are defined (analogous to light) as the states with 
a variance in one quadrature exciton operator smaller than that associated with the 
coherent state, the exciton becomes squeezed if 

<:(AX(t)) 2 :> = < (AX(t)) 2 > -\ < 0, (15) 

or <: (AP(t)) 2 :> = < (AP(t)) 2 > ~ < 0, (16) 

where <: (AX(t)) 2 :> and <: (AP(t)) 2 :> are so-called the normally ordered variances 
and are widely accepted as a measure of squeezing [15,16]. 

Let us, for simplicity, suppose that at initial time t — o excitons are in coherent state 
and light is prepared in a squeezed state characterized by the squeezed factor r and a 
phase if [1]. From (9), (12)-(14) and after some algebraic manipulations we obtain the 
following exact analytical expressions for the exciton normally ordered variances 

<: (AX(t)) 2 :> = \\(3'{t) | 2 sinh 2 r-^Re(/3 ,2 (t)e^) coshrsinhr, (17) 
<: (AP(t)) 2 :> = \\(3\t) | 2 sinh 2 r + ^Re(/3' 2 (t)e^) coshrsinhr, (18) 
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where 

|^'(t)| 2 = 4 M Ksin 2 ^^, (19) 

Re{(3' 2 {t)e 1 ^) = -Au\v\ sin 2 (ei - ^ cos(e 1 + e 2 - <p)t. (20) 

where ei,e 2 and Ui,vi are defined above by (6)-(8) and (10), (11), respectively. 

The formulas (17)- (20) are our main results from which we would like to emphasize 
the followings. Since u i: Vi are expressed in e i: E* and g* (see (10), (11)) which in turn 
depend explicitly on light frequency ui , coupling constants /, A , and the exciton density 
n ex (see (6)- (8)), one can expect various kinds of dependence of the exciton squeezing 
on these parameters. 

Note that the exciton density n ex always comes along with the parameters / and 7 
characterizing the nonlinear exciton-exciton interaction and anharmonicity of exciton- 
photon coupling due to phase space filling of exciton. The coupling constants / and 7 
were calculated for the lowest lS-exciton in bulk semiconductors in [9,14] 

f = ^E B a%, 7 =^4. (21) 

where a B is the effective Bohr radius of the bulk exciton and E B is its binding energy. 

For graphical illustration we take the case of CdS with the parameters E ex = 
2553mev , g = 50mev , a B = 2.8 x 10~ 7 cm , E B = 27.8mev and (p = 0, r = 0.8, 
scaled time is equal to E ex t. To see the frequency dependence of the exciton squeezing, 
in Fig. la, b, c we plot the time evolution of the exciton normally ordered variance 
<: (Ax(t)) 2 :> for the fixed exciton density n ex = 1.10 17 cm -3 , and three different 
values of photon energy Auj = E ex — u = 350meV, 100 meV and meV, respectively. 
One can see that with photon frequency approaching closer to the exciton resonance 
( Auj — > 0) the squeezing degree obtained by excitons from photons effectively increases, 
and gains its maximum at exact exciton resonance (Fig. lc). 

One of the new main results in this paper is explicit dependence of exciton squeezing 
on the exciton density. This is a consequence of the effects of nonlinear interaction 
between excitons and the anharmonic exciton-photon interaction due to phase space 
filling of excitons. Higher exciton densities more screen the interaction between excitons 
and photons, and reduce it by the factor (1 — n ex VX). This is just the saturation effect 
due to phase space filling of excitons which usually appears in highly excited systems. 
The above-analyzed density dependence of exciton squeezing is well seen from Fig. 2a, 
b, c plotted for the fixed photon frequency near exciton resonance Auj = 3 meV, and 
some different exciton densities n er . 
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Figure Captions 



Fig.l Time evolution of the exciton normally ordered variance <: (Ax(t)) 2 :> for 
n ex = 1.10 17 cm' 3 ; Au = 350meV (a), 100 meV (b) and meV (c). The closer to 
exciton resonance the greater exciton squeezing achieved. 



Fig. 2 Time evolution of the exciton normally ordered variance <: (Ax(t)) 2 :> for 
Alo = 3meV; n ex = 1.10 17 cm' 3 (a), 1.10 18 cm' 3 (b) and 3.10 18 cm" 3 (c). With 
increasing n ex the exciton squeezing degree decreases. 
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